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Abstract. We show how two-point correlation functions derived within non- 
isotropic random wave models are in fact quantum results that are obtained in 
the appropriate limit in terms of the exact Green function of the quantum system. 
Since no statistical model is required for this derivation, this shows that taking the 
wave functions as Gaussian processes is the only assumption of those random wave 
models. We also show how for clean systems the two-point correlation function 
defined through an energy average defines a Gaussian theory which substantially 
reduces the spurious contributions coming from the normalisation problem. 

PACS numbers: 05.45.Mt,05.40-a 

Since Berry's seminal paper in 1977 the so-called Random Wave Model 
(RWM) has become by far the most popular and successful tool to describe the 
statistical properties of wave functions of classically chaotic systems which in this 
approach are modelled by a random superposition of plane waves. Its applications 
range from the realm of optics [2], passing by the general problem of wave mechanics 
in disordered media |3| to important issues in mesoscopic systems 0]. Owing to this 
robustness this approach has been regarded as the indicator of wave signatures of 
classical chaotic dynamics f^- 

The reasons for this success can be traced back to two fundamental points. First, 
it can be formally shown that such a random wave function is a stationary random 
process |S] (roughly speaking a function taking random values at each point); second, 
such random process is Gaussian, namely, it is uniquely characterised by a two-point 
correlation function which expresses fundamental symmetries, like the isotropy of free 
space. The fact that the process is Gaussian represents a considerable advantage in 
an operational sense since it provides us with a set of rules to cope with averages 
over complicated expressions in the way Wick's theorem and its variants do. At the 
same time the generality of the random wave two-point correlation makes the theory 
a remarkably good approximation when the effect of the boundaries can be neglected, 
like in the case of bulk properties. 

When applied to real quantum systems, however, there remain limitations related 
to the above mentioned ingredients. Goncerning the Gaussian assumption, a formal 
proof showing that a chaotic wavcfunction is indeed a Gaussian process is still 
lacking. Even more, as noted in 7 the Gaussian distribution explicitly contradicts 
the normalisation condition for the wavefunction. In practical terms this means 
that, when dealing with statistics beyond the two-point correlation function, the 
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Gaussian distribution produces spurious non physical contributions, and attempts to 
construct a RWM respecting the normaUsation constraint lead to severe mathematical 
difhculties |S] . Still this assumption is supported by many arguments based on Random 
Matrix Theory 0] , quantum ergodicity |S] , information theory 8 , and Berry's original 
semiclassical picture Impressive numerical results also support the conjecture at 
the level of one-point statistics [HI, and evidence for higher order statistics is given 
in 1^12]. Hence, it is appealing to look for a RWM which minimises the spurious 
contributions due to the normalisation problem while keeping the wavefunction 
distribution still Gaussian. 

About the isotropic character of the theory, constructing a random superposition 
of waves satisfying both the Schrodinger equation and boundary conditions turns out 
to be at least as difficult as solving the full quantum mechanical problem by means of 
standard techniques. To our knowledge the attempts in the direction of a non- isotropic 
RWM can only deal with highly idealised boundaries such as an infinite straight wall 
|1U|. a linear potential barrier (11) . and the edge between two infinite lines enclosing an 
angle of a rational multiple of tt • Also in ^] a variation of the RWM to include 
finite size effects is presented. The fact that these approximations already produce 
non-trivial deviations form the isotropic case is an indication of the importance of the 
inclusion of arbitrary boundaries. 

Our aim in this communication is twofold: first, we shall show that the mentioned 
results for the two-point correlation function defining the non-isotropic and finite- 
size RWM can be derived from quantum mechanical expressions, namely, they are 
independent of any statistical assumption about the wavefunction. Second, we shall 
show how for a statistical description of wavefunctions using an energy ensemble 
average, the spurious contributions coming form the normalisation problem are of 
order 0{1/N) with N the number of members of the ensemble, making their effect 
negligible for high energies. 

Isotropic and non-isotropic random wave models. We consider solutions of the 

Schrodinger equation + V{r)j '0„(r) = Eni^ni"^) for a closed system where 

the corresponding classical dynamics is chaotic (in the following we take 2m =1). The 
RWM assumes the statistical description of an ensemble of wavefunctions mimicked by 
a random superposition of plane waves with local wavenumber k{r) = \/e — V{r)/h, 
where e is the mean energy of the states under study. For the sake of comparison we 
shall focus on the following averages, used for the nodal counting statistics (we follow 
the notation of Berry 10 ): 



We wish to stress, however, that the RWM, being Gaussian, can deal with far more 
general averages. For billiard systems the isotropic RWM (denoted by a superscript 
i) is defined by the ensemble 



and the average (. . .) is defined by integration over a set of independent random phases 




(1) 
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e (0, 2tt]. We also choose 9j — 2Trj/J, where the hmit J oo is always taken after 
averaging over the 0's. Explicit calculation then gives the following results 

i?^(f) = 1, Dlif) = ^, Dlif) = ^, Klir) = 0. (3) 

These results represent bulk approximations to the real situation, since boundary 
effects are completely neglected. In order to improve this limitation, the following 
ensemble of non-isotropic (ni) superpositions of random waves was introduced in lO] 
to take into account the effect of a straight infinite boundary aX y — on which 
we demand the wavcfunction to satisfy Dirichlct (D) or Neumann (N) boundary 
conditions: 

— ^ sin k{y — yo) cos{kx cos 9j + (l>j), (4) 

ip^{r) = \ -f ^ cosk{y — yo) cos(fcxcos0j + (l>j). (5) 

With the averaging procedure as in the isotropic case, one obtains for the Dirichlet 
(upper sign) and Neumann (lower sign) cases [111] : 

B"\r) = It Jomy-yo)), 

DTir) = y(1 T Jo(2fc(y - yo)) T J2(2A:(y - yo))), (6) 

A7« = y(1 ± Jomy - yo)) T J2{2k{y - yo))), 
K^\r)^±kJ,{2k{y-yo)). 

For more general situations where the confining potential is smooth (S), Bies and 
Heller introduced the following ensemble of random Airy functions Ai(x) to satisfy 
locally the Schrodinger equation for a linear ramp potential V{x, y) = Vy: 

1 ' 

^Hr) ^^J2Ai [^{y, Q,)] exp [i{Q,x + cb,)]. (7) 

Here 



^{y,Q)={Pl\y~yo)+(yy Q', 

and yo = e/V is the turning point in the direction of the linear ramp, fixed by the 
mean energy e of the eigenstates under study. The phases (jjj are defined as usual and 
provide the averaging, while Qj £ [— oo,cx)]. Explicit calculation then gives | 11[I14| : 



B^(r) = 


/ Ai' [^{y,Q)]dQ 

JO 




Dlir^ = 


/•OO 

/ [*(?/, Q)]dQ 

Jo 


(8) 


D'yif) = 


/ Ai'^ [^{y,Q)]dQ 
Jo 




K^ir) = 


/ Ai[vI/(y,g)]Ai'[vI/(y,g)]dQ 
Jo 
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where Ai'(a;) is the derivative of the Airy function. 

The quantum description. We consider a set of normahsed solutions ipn{r) of 
the Schrodinger equation with non-degenerate eigenvalues lying inside the interval 
W = [e — ^,e+ Considering W as a range of energies with almost constant mean 
level spacing A(e), the number of states within the interval is iV = (in general 

= si p{E)dE with p{E) the level density). In the high-energy limit we are 

interested in, ^ 1 with ^ <C 1 are well defined limits which we shall always assume 
implicitly. The two point correlation function, 

F{r^ur2):^mriW{r2)):=^ E Mri)rn{r2) 

can be used to calculate the averages in Eq. (1) by differentiation: 

[F{ri,r2)]p^^P^^p, 



(9) 



B{r) := 
l?.(r) 



Ky{r) 



dxidx2 
92 



F{ri,r2) 
F{fi,r2) 



(10) 



Firi,r2) 



_dyidy2 
"1 _d_ 

2 \dyi dy2. 
It is convenient to use the Green function of the system, 

i^n{riWn{r2) 



Giri,f2,E + iO+) = 



^ E- 

n=l 



En + iO+ ' 



(11) 



to obtain the expression 

F{ri,r2) 

A(e) 1 
27ri Se 



{G* {n,r2,E + iO+) - G(f2 , fi , ^ + iO+))dE. (12) 



Note that this is an exact result and the common approximation F {7^1,7^2) 
{G* {fi,f2,e + iO"*") — G{f2 J ) e -|- iO"*")) |15| is not valid in general and requires further 
assumptions. Even more, the additional energy integration in Eq. (12) will turn out 
to be essential. 

Different approximations to the Green function valid under different situations 
can now be used to study the corresponding wavefunction statistics. 

The bulk contribution and finite size effects. For billiard systems the bulk (6) 
results are obtained by using the free propagator given in two dimensions by the 
Hankel function G"{f'2,ri, E + iO+) — ^77g^^ (-^jri — ^2 1) instead of the exact Green 
function. The corresponding contribution to the two-point correlation is: 



\ri - r2\ 



dE 



(13) 



where A is the billiard area. Using Eqs. (10) and (13) we easily recover the results in 
Eq. (3) correctly normalised. Further analysis of expression (13) shows that it reduces 
to Berry's result P 



r2 



(14) 



Letter to the Editor 



5 



when \ri — ^ \/4:A/tt, while it decays much faster for jri — ^2! > y^4A/TT as long 
as Se > hy/ire/AA. Noticing that y/AA/n is just the average system linear size L we 
see that Eq. (13) actually defines a RWM which incorporates finite size efl^ects when 
the average is taken on scales larger than the ballistic Thouless energy, eth = h^/e/ L. 
Eq. (13) then provides an analytical expression for the correlation function defined in 

The case of an infinite straight barrier. For this situation we construct Green 
functions with the correct parity under the reflection symmetry with respect to the 
line y ~ yo hy means of the method of images. The symmetric and antisymmetric 
combinations give the corresponding two-point correlation for Dirichlet (upper sign) 
and Neumann (lower sign) boundary conditions as: 



(^1,^2) 



1 1 

'aJ^ 



\/{xi - X2Y + (yi - 2/2)^ 



E 



X2Y + (2yo - 2/1 - yif 



dE. 



(15) 



Using this correlation function and Eq. (10) we obtain the averages defined in Eq. 
(1). Berry's results Eq. (6) are again obtained in the limit of very short distances to 

the boundary |y — yol ^ 



The infinite, smooth barrier. For a particle in the potential V{x, y) = Vy, the 
Schrodinger equation is separable. The solutions along the x direction are plane waves 
and in the y direction Airy functions. Hence we have: 



ikkAx, y)= \ I exp {~ikx)M 



(-7) 



and for the Green function 

G(fi,f2,S + iO+) 



00 />oo 



— 00 <J —00 



E-e- h?B + ^0^ 



-dedk. 



Including this result into Eq. (12) we find for the correlation function: 



cos {k{xi — X2)) 



Ai 



^) 12/1 



E - h^k"' 
V 



Ai 



E - h^e 

V 



(16) 
(17) 

(18) 
dkdE. 



Together with the relations Eq. (10) again, one obtains Eq. (8) in the limit A(e) 0. 

To summarise so far, we obtained the one-point averages, Eq.(l), for a closed 
system from pure quantum mechanical considerations without appealing to any 
statistical assumption about the wavefunction. The known RWM results Eqs. (3,6,8) 
are then derived in the appropriate limits (either short distances or infinite system 
size) . The generalisation to any other average bilinear in the wavefunction is straight 
forward. 

The normalisation problem. There is a prominent counter argument against the 
Gaussian assumption first presented to our knowledge in |3| and further explored in 
which deserves special attention. It is the apparent contradiction between the 
normalisation of the members of the ensemble and the Gaussian distribution of the 
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wavefunction amplitudes. Mathematically this can be stated in the following way. 
Consider the functional 

m = J mrll^df, (19) 

where ip{r) is a member of the ensemble we use to describe the statistical properties 
of the wavefunction. One constraint we must impose is the normalisation of all 
expressed by the vanishing of the ensemble variance Var{ri) — {{r][tp])'^) — (('7[V-']))^- 
If the wavefunction's distribution is Gaussian, we obtain 

Var{r,)^2[ ( \{WiW {r2))?dridr2- (20) 



This is in clear contradiction to the normalisation condition Var{rii) = 0. Here we show 
that this variance is of order Var{ri) = 0{1/N). To this end we recall the definition 
of the two-point correlation, Eq. (9), and note that the ipnir) are eigenfunctions of 
the same Hamilton operator, i.e, they form an orthonormal set: 

V'^(r)^*(r)(if-(5,,,. (21) 

Then it is easy to obtain the following composition rule for the two-point correlation: 

{^P{n)r{r^){Hr^r{r2))dr = l(^(ri)^*(f2)). (22) 

Since {tp{r{)^* {f2)) converges in the Umit ^ ^ const, ^ oo, we see that indeed 
Var{rj) converges to zero as 0{\/N). 

This behaviour of Var{ri), relying on the fact that an energy ensemble average is 
taken, differs from the case of disordered systems where the lack of orthogonality 
between the different members of the ensemble (since they are eigenfunctions of 
different Hamiltonians corresponding to different disorder realizations) leads to an 
expression for Varirf) of order 0(1). 

It is important to note that the extra energy average is essential to satisfy Eq. 
(23). For example, a simple calculation shows that Berry's result, Eq. (14), does not 
satisfy the composition rule, while our result Eq. (13) does it as long as 5e > exh- This 
is a particular case of a more general statement saying that under certain conditions 
approximate Green functions will make the correlation function satisfy the composition 
rule. The proof of this result requires the use of semiclassical techniques and will be 
presented elsewhere [T7] . 

Concluding remarks. We have shown that all the two-point correlation functions 
used to fix the different random wave models, (isotropic, non-isotropic, for a smooth 
boundary and including finite size effects) can be derived in the appropriate limit of 
the exact quantum mechanical expressions. To this end we do not use any statistical 
assumption about the wavef unctions, in fact, this results are independent of the regular, 
mixed or chaotic character of the classical system |1()| . Also, we showed that for clean 
chaotic systems the use of an energy ensemble reduces the spurious contributions 
coming from the normalisation problem without affecting the Gaussian assumption. 
This result only requires the consistent use of the quantum mechanical definition of 
the correlation function. 

A Gaussian field with a correlation given in terms of the exact Green function 
of the system and the energy average carefully taken into account represents a 
generalisation which includes all known RWM as limiting cases and successfully takes 
into account boundary and normalisation effects for any closed, clean chaotic system. 
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To deal with general shapes (like the ones study in ^2]) and boundary conditions (like 
the mixed case presented in ^3]) , however, the exact quantum approach presented 
here cannot be analytically performed. An approach using the semiclassical Green 
function to derive the asymptotic expressions for the correlations presented here (and 
the ones in ^1 is the adequate method and will be the subject of a separate 
communication |17j . 
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